Non-Hamiltonian Commutators in Quantum Mechanics 
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The symplectic structure of quantum commutators is first unveiled and then exploited to introduce 
generalized non-Hamiltonian brackets in quantum mechanics. It is easily recognized that quantum- 
classical systems are described by a particular realization of such a bracket. In light of previous work, 
this introduces a unified approach to classical and quantum-classical non-Hamiltonian dynamics. In 
order to illustrate the use of non-Hamiltonian commutators, it is shown how to define thermodynamic 
constraints in quantum-classical systems. In particular, quantum-classical Nose-Hoover equations of 
motion and the associated stationary density matrix are derived. The non-Hamiltonian commutators 
for both Nose-Hoover chains and Nose- Andersen (constant-pressure constant temperature) dynamics 
are also given. Perspectives of the formalism are discussed. 
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I. INTRODUCTION 

In order to describe phenomena in the real world, clas- 
sical and quantum systems are represented by means of 
Hamiltonian mathematical theories P, However, 
when studying systems with many degrees of freedom, 
the need of performing numerical calculations on com- 
puters has led to the development of non-Hamiltonian 
mathematical structures |5|, L* Li • 

In the classical case, non-Hamiltonian formalisms are 
typically employed to implement thermodynamic con- 
straints [8j using just few additional degrees of freedom 
(whereas by using Hamiltonian theories one should re- 
sort to an infinite number of degrees of freedom). Just 
recently, it has been shown that the non-Hamiltonian 
dynamics of classical systems can be formulated in a uni- 
fied way by means of generalized brackets which ensure 
energy conservation [£| E3 and subsumes Dirac's for- 
malism 3] for systems with holonomic constraints 
Other approaches to classical non-Hamiltonian brackets 
can be found in Refs. [T^. [l3| . 

In the quantum case, the impossibility to solve on 
computers full quantum dynamics for interacting many- 
body systems has led to the development of quantum- 
classical theories. Indeed, a generalized bracket to treat 
quantum-classical systems has been proposed by various 
authors fl4j. 

Since in the classical case non-Hamiltonian brackets 
are obtained by modifying the symplectic structure of the 
Poisson bracket 0, 0] , in order to deal with the quan- 
tum case one could first make apparent the symplectic 
structure of the commutator (which is the Hamiltonian 
bracket of quantum mechanics) and then generalize it 
in order to obtain a non-Hamiltonian quantum bracket 
(commutator). In this paper it is shown that this is in- 
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deed possible. The non-Hamiltonian commutator, which 
is obtained by this procedure, is then used to reformu- 
late quantum-classical brackets Thus, it is stressed 
that quantum-classical dynamics can be regarded as a 
form of non-Hamiltonian quantum mechanics because the 
quantum-classical bracket does not satisfy the Jacobi re- 
lation and, as a consequence, the time-translation invari- 
ance of the algebra is violated. In order to illustrate the 
use of non-Hamiltonian commutators, it is shown how to 
define thermodynamic constraints in quantum-classical 
systems. The particular case of the Nose-Hover thermo- 
stat 0, is treated in full details and the associated sta- 
tionary density matrix is derived. The more general cases 
of Nose-Hover chains 0] and constant pressure and tem- 
perature bring no major difference neither con- 
ceptually nor technically and are treated in less detail. It 
is worth noting that some past attempts of introdu cing 
Nose-Hover dynamics in quantum calculations [Tfl Il7| 
used a simpler form of quantum-classical dynamics which 
did not treat correctly the quantum back-reaction on 
the classical variables. The possibility of applying ther- 
modynamic constraints to quantum-classical dynamics is 
a technical advance that could lead to further theoreti- 
cal and computational achievements with regards to the 
study of open quantum systems 0, 0] . With respect 
to this, a non-trivial major obstacle is the development 
of efficient algorithms to simulate long-time quantum dy- 
namics. 

Besides the technical applications of non-Hamiltonian 
commutators to the particular case of quantum-classical 
dynamics, one could appreciate on a more conceptual 
level that, in light of previous work, non-Hamiltonian 
brackets provide a unified approach to non-Hamiltonian 
dynamics both in the classical and quantum case. In 
addition, if one is willing to indulge in speculations, it is 
worth to note that the mathematical structure presented 
in this paper may be shown to generalize the formalisms 
that a number of authors have already presented in the 
literature [U H III HI HI In particular it is worth 
mentioning that non-Hamiltonian commutators could be 
used, in principle, in order to introduce non-linear effects 
in quantum mechanics along the lines already proposed 
by Weinberg j^]. Therefore, one could foresee interesting 
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applications of non-Hamiltonian commutators in various 
fields. 

The paper is organized as follows: in section [H] the 
symplectic structure of Hamiltonian quantum mechanics 
is unveiled and its generalization by means of the non- 
Hamiltonian commutator is proposed. In section II I II it 
is shown that the quantum-classical bracket can be writ- 
ten in matrix form as a non-Hamiltonian commutator. 
Such a form easily illustrates the failure of the Jacobi 
relation [25| . In section IIVI non-Hamiltonian commu- 
tators for quantum-classical systems are used in order 
to introduce, following Refs. 0, EJ, Nose thermostat- 
ted dynamics on the classical degrees of freedom. Nose- 
Hoover chains and constant pressure constant tempera- 
ture (NPT) equations of motion bring no major difference 
and are treated in less detail in appendix lAl In scctionlVl 
it is proven that the stationary density matrix under the 
quantum-classical Nose-Hoover equation of motion exists 
and its explicit form is given up to order h. Conclusions 
and perspectives are given in the final section. 



II. NON-HAMILTONIAN QUANTUM 
MECHANICS 



It is well-known that classical and quantum dynam- 
ics share an analogous algebraic structure 0, Q realized 
by means of specific brackets: Poisson brackets in the 
classical case and commutators in the quantum one. It 
is also known that Poisson brackets have a symplectic 
structure that is easily represented in matrix form Q, • 
If one denotes the point in phase space as X = (R,P), 
where R and P are the usual coordinates and momenta 
respectively, by defining the antisymmetric matrix 



B 



1 

-1 



(1) 



the Poisson bracket of any two phase space function a(X) 
and b(X) can be written in matrix form as 



2N 



da 



i)h 



M= L QX^dX 



(2) 



where 2N is phase space dimension. In Refs. |^. ITollrT| . 
Eq. J5J has been generalized introducing an antisym- 
metric tensor field Bij(X) = —Bji(X) so that a general 
bracket {...,.. .}x, having the same matrix structure of 
that in Eq. J2J), could be introduced and non-Hamiltonian 
equations of motion could be written as 



2iV 



X i = {X i ,H}x = Y,Bi j (X) 

3 = 1 



an 

8X~ 



(3) 



where Ti is the "Hamiltonian" or generalized energy. 

As one could expect, the commutator can also be 
written in matrix form using the symplectic structure 



of Eq. Q). If one considers a set of quantum vari- 
ables x a , a — l,...,n, which can be canonical, non- 
canonical or anti-commuting variables, the commutator 

[Xa,Xu] = XaXu~XuXa (a, v = can be expressed 

as 



[Xa,X»] = [Xa Xu] 



o r 




Xa 


-1 







(4) 



The above matrix form of the commutator permits to ap- 
preciate the common symplectic structure of both classi- 
cal and quantum mechanics. 

Given the Hamiltonian operator H of the system, the 
law of motion in the Heisenberg picture can also be writ- 
ten in matrix form as 



dXa 
dt 



[ H x« ] ■ B ■ 



H 

X a 



= iC-Xc 



where it has been introduced the Liouville operator 



iC 



h 



[H 



B 



H 



(5) 



(G) 



The algebra of commutators is a Lie algebra. This 
means in particular that the commutator satisfies the 
following properties: 

\Xa,Xv] = -[Xv,Xa\ (7) 
[XaXviXo] = Xa[Xu,Xa] + [x a ,Xa]Xv (8) 

[c,Xu] - 0, (9) 

where c is a so called c-number and a,v,a = 1, ...,n. 
Besides properties in Eqs. (|7I9|I . in order to have a Lie 
algebra, it is necessary that the so called Jacobi identity 
holds 

J = [Xa, [Xu, Xa]] + [Xv, [Xa, Xu]] + [&/, [tier, Xa]] = 0. 

(10) 

The Jacobi identity ensures that the algebra is invariant 
under the law of motion and as such it states an inte- 
grability condition. In the above formalism it can be ap- 
preciated that the antisymmetry of the commutator Q 
arises from the antisymmetry of the symplectic matrix 
B and ensures that if H is not explicitly time-dependent 
then it is a constant of motion 



d 
dt 



H = iCH = 



(11) 



The conservation of energy under time-translation de- 
fined by means of antisymmetric brackets is another nice 
property shared both by the algebra of Poisson brackets 
on classical phase space and by the algebra of commuta- 
tors of quantum variables. 

Using the operator language of Eq. one can define 
a generalized commutator as 



\Xa,Xv] = [Xa Xu]-T> 



Xa 
Xu 



(12) 
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where 'D is an antisymmetric matrix operator of the form 



o C 
-C o 



(13) 



with £ arbitrary operator or c-number. Generalized equa- 
tions of motion could then be defined as 



dt 



k 



[ H Xa ] • T> 



H 

Xa 



it Xa - (14) 



It must be stressed that the non-Hamiltonian commu- 
tator defined in Eq. (|12fl could violate the Jacobi re- 
lation H1(J|) so that in general it does not define a Lie 
algebra. The non-Hamiltonian commutator of Eq. 112|) 
defines, of course, a generalized form of quantum mechan- 
ics. However, in this generalized theory the Hamiltonian 
operator H is still a constant of motion because of the an- 
tisymmetry of T>. It is interesting to note that T) could in 
principle depend from the quantum variables Xa- Then 
Eq. (|14fl can be thought of as a generalization to the 
Heisenberg picture of the mathematical formalism pro- 
posed by Weinberg p3| in order to introduce non-linear 
effects in quantum mechanics. 

In the next section it will be shown that the non- 
Hamiltonian commutator defined in Eq. 112|) and the 
non-Hamiltonian equations of motion l|14|) provide the 
mathematical structure for quantum-classical evolu- 
tion pn. 



Exploiting what has been done in Refs. 0, for the 
Poisson bracket and in the previous section for the com- 
mutator, the quantum-classical bracket can be easily re- 
casted in matrix form as a non-Hamiltonian commutator. 
To this end, one can introduce the operator A defined 
in such a way that applying its negative on any pair of 
quantum-classical operators functions X a{X) and x v [X) 
their Poisson bracket is obtained 



27V „„ 

{Xa,X,} = -X a (X)A X ,AX) = 



d X u 
dXi 



(16) 



The quantum-classical law of motion can be rewritten 

as 



dtXa = j: [ H Xa ] ■ B ■ 



[H 



H 

Xa 



r o - ' 

A 2 




' H ' 


- A 
L 2 u _ 




Xa 



(17) 

A more compact form is readily found by defining the 
antisymmetric matrix super-operator 





(1 + ^) 



1 



HA 
2 1 







(18) 



Using the matrix super-operator in Eq. 1)18(1 the 
quantum-classical law of motion becomes 



III. NON-HAMILTONIAN COMMUTATORS IN 
QUANTUM-CLASSICAL MECHANICS 

Quantum-classical systems can be treated by means of 
an algebraic approach. This has been already proposed 
by a number of authors |14| by means of a quantum- 
classical bracket which does not satisfy the Jacobi rela- 
tion. A quantum-classical system is composed of both 
quantum X an d classical X degrees of freedom. The 
quantum variables depends from the classical point X 
so that an abstract space is defined in such a way that a 
Hilbert space (where quantum dynamics takes place) is 
attached to each phase space point. In turn, a displace- 
ment of the phase space point determines a consistent 
effect on quantum evolution in the Hilbert space. The 
energy of the system is defined in terms of a quantum- 
classical Hamiltonian operator H = H(X) coupling 
quantum and classical variables E — Tr' J dXH(X). It 
has been shown [l4| that the dynamical evolution of a 
quantum-classical operator x(A) is given by 



dtXa = J. [ H Xa]-V- 



H 

Xa 



= (H,Xa)=*£Xa, (19) 



d tX (X) = -[H,x(X)} 
= (H,X(X)). 



l -{H,x{X)}+ l -{x{X),H} 
(15) 



The last equality defines the quantum-classical bracket 
in terms of the commutator and the symmetrized sum of 
Poisson brackets. 



where the last equality introduces the quantum-classical 
Liouville operator in terms of the quantum-classical 
bracket. The structure of Eq. (|19[) is that of the non- 
Hamiltonian commutator given in Eq. I|14[) and as such 
generalizes the standard quantum laws of motion of 
Eq |J5J|. It is clear from its definition in Eq. (|18fl that the 
antisymmetric matrix super-operator T> has not a sim- 
ple symplectic structure as 13. It contains the operator 
A defined in Eq. H16fl which, in this case, has a sym- 
plectic structure. As such T) introduces a novel math- 
ematical structure that characterizes the time evolution 
of quantum-classical systems. 

The Jacobi relation in quantum-classical dynamics is 

J = (Xa, (Xv,X<j)) + (X<j, (Xa>Xu)) + (Xv, (X<r,Xa))- (20) 

Using the matrix formalism introduced it is simple to 
calculate J explicitly and to this aim one can consider the 
first term on the right hand side of Eq. H2U|) . The other 
two terms on the right hand side of Eq. I|20() can then be 
easily calculated by considering the even permutations 
of Xa, Xu, Xa m the formula obtained for ( X a, {XviXa))- 
Finally, collecting the terms together one gets 



J = -g{Xak{Xukx<?) - XaA(x CT Ax„) - {Xu^X<y)J^Xa 
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+ (x<x A X^) A Xa + X CT A(x Q Ax„) - XuHXa^Xa) 
~ (Xa^Xv)AXo- + {Xa^X^^Xu + XvMXvA-Xa) 
~ XaMXv^Xa) ~ (XakXa)AXv + (X^ A Xa) A X<x} • 

(21) 

In order to easily get such expression the following rela- 
tion 

2N 2N 



dXc 
dXi 



+ ig-BijMOiX,) (22) 



was exploited. Thus it is found that the Jacobi relation 
does not hold globally for all points X of phase space 
(J * 0). 



IV. NOSE DYNAMICS IN 
QUANTUM-CLASSICAL SYSTEMS 



The antisymmetric matrix 3 enters through A in the 
definition of T>. Following the work of Ref. 0,^3 thermo- 
dynamic constraints can be imposed on the classical bath 
degrees of freedom in quantum-classical dynamics just 
by modifying the matrix 3. For clarity it will be explic- 
itly shown how to generalize the derivation of quantum- 
classical equations of motion in the case of Nose constant 
temperature dynamics . In the Nose case the bath de- 
grees of freedom will be 



X=(R,r),P, Pri ) , 



(23) 



rj and p v are the Nose coordinate and momentum. The 
following quantum-classical Hamiltonian is assumed 



The explicit form of the matrix operator that defines 
through Eq. (|19(1 the quantum-classical bracket and the 
law of motion is then given by 







l + ^A N 

6 



(27) 



by 



The quantum-classical Nose-Liouville operator is given 



X 



(28) 



One is then led to consider, in the right hand side of l|28|l . 
the term given by 



- HnAnX + xA N i? N = 



- 2F, 



<9$ dx d X <9$ 
dRdP + &P~dR 



-2 — — - o p v dx 
M dR m v dr\ 

+2^P§ , 
m v aP 



dx 

V dPr, 



(29) 



p 2 

where F v = jg- — gkgT . Finally using the above result 
the equation of motions for the dynamical variables are 
given by 



dt 



X 



h \ x x > 2 \dPdR dRdP 



P dx , Pn dx p v p dx 



F, 



dx 



M dR m v di] m v dP v dprj 



(30) 



A. Representation in the Adiabatic Basis 



K 



p2 



P 



2M 2m, 



■*(jc,R) + gk B Tri, (24) 



where K is the quantum kinetic operator, $ is the poten- 
tial operator coupling classical and quantum variables, 
M is the mass of the classical degrees of freedom, m n is 
Nose inertial parameter, and g is a numerical constant 
whose value (as it will be shown) must be set equal to 
the number N of classical momenta P if one wants to ob- 
tain a sampling of the R, P coordinates in the canonical 
ensemble. Then the matrix B N is 



3 



10 

1 

-10 -P 

-1 P 



(25) 



Using 3 the operator An and the classical phase space 
non-Hamiltonian bracket on two generic variables A\ and 
A 2 can be defined 



A X A N A 2 



2N 



dA 1 N dA 2 



(26) 



One can express the quantum-classical equations of 
motion in the adiabatic states. Nose quantum-classical 
Hamiltonian can be rewritten as 



h(R) 



P 2 , V% 



2M 2m 



gk B Tr] , (31) 



where it has been introduced the operator h(R) = K + 
$(x, R). Then the adiabatic states are defined by 

h(R)\a;R) = E a (R)\a; R) . (32) 

In the adiabatic states, Eq. (fTTTTf) is easily found to be 



d 1 1 P 8v aa 

u aa ■ aa , 

dt X = ^' X + M^R 



-P 



Pv 8 p v d 



~^+Fn 



d 



M 

2 dP 



m v dP rriri drj v dp n 
P 

Ti 



x 



P . 0a' P a/3' j , 1 dx"" 13 T^B'a' 



2 dP 



(33) 
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where F a ^ 



9* I 
077 I 



and g? q/ 3 = (a\-£>\p) is the 
nonadiabatic coupling vector. Equation (|33l) can be 
rewritten introducing the Liouville operator iC N such 
that 



E 



13/3' 



(34) 



The operator is 



,-pN 



— iuj aa '5 a p5 a > pi + 5 a p5 a i(j 
\ m v dP 



P d 
JldR 
Pn d 

-5 aS F p ' a ' 'A + -5 a ,B>F af} -^- + —d ad 5 a ,g> 
2 a P dP 2 p dP M p p 

P 

jjd/3' a i5 a p . (35) 



The quantum-classical Liouville operator can be put into 
a form that makes its structure more apparent by adding 
and subtracting the term 



V. STATIONARY NOSE DENSITY MATRIX 



The average of any operator \ can be calculated from 

(X) = Tr' J dX p N x(t) = Tr' / dX p N exp (iC N t) X . 

(41) 

The action of exp («£ N i) can be transferred from \ to fa 
by using the cyclic invariance of the trace and integrating 
by parts the terms coming from the classical brackets. 
One can write 



~\{ 11 {..-.//n}). (42) 



In this equation the classical bracket terms are written 



2JY 



]//x....i !....//x! = E i^Bif 9 ' 



5 a p5 a 



Then using 



F 



a/3 



F a + (E a — Ep) d a p 



(36) 



(37) 



and rearranging the terms one obtains a classical-like 
Nose-Liouville operator 



P d 
MdR 



pVrt_d_ + Pr, 9 



+ ^F a + F a ' 
and a jump operator 

— Jota',/3/3' = fia'/3'd a p 
+ Safid^p, 



rrir! dP 
d_ 
dP 



F, 



8 



rrir, drj 11 dp v 



(38) 



P 

M 
P 
M 



E a -E fi } d_ 
dP 



in terms of which the quantum-classical Liouville opera- 
tor is finally written as 



lJ -aa',p0' 



lLO aa i5 a p8 a 'p' + S a p5 a 'p'iL aa , — Jaa',/3/3'- 

(40) 

The jump operator Jaa',/3p' is responsible for transitions 
between adiabatic states while the classical-like Nose Li- 
ouville operator iC^ a , pgi expresses Nose dynamics on 
a constant generalized energy surface with Hellman- 

Feynman forces given by 1/2 i^F a + F a J . 

This shows that the matrix form of the non- 
Hamiltonian commutator is suitable for the development 
of generalized non-Hamiltonian dynamics for classical de- 
grees of freedom in quantum-classical systems. 



dXi » dXi 



d... N dH N 
dX t 13 dXi 



(43) 



When integrating by parts the right hand side, one ob- 
tains a term proportional to the compressibility kn = 

Si 7=i TJx^" ^ s a resmt the m i xe d quantum- 

classical Liouville operator, in this case, is not hermitian 

(iC N ) f = -tC N - k n • (44) 

The average value can then be written as 

(X) = Tr' / dX xexp [-(iC N + « N )t] p N (45) 



The mixed quantum-classical Nose density matrix 
evolves under the equation 



d 



i 



H N , PN 



±({H N ,fa}-{fa,H"} 



(46) 



The stationary density matrix fa c is defined by 

(i£ N + KN )p Ne - . (47) 

To find the explicit expression one can follow Ref. |25| . 
expand the density matrix in powers of Ti 



71=0 



(48) 



and look for an explicit solution in the adiabatic basis. 
In such a basis the Nose-Liouville operator is expressed 
by Eq. (|40|l and the Nose Hamiltonian is given by 



Hi 



p2 . .'-> 

2M ' 2m 



^T7 + ^+9kBT V + E a (R) 



K(R,P) 



9 
Pt, 

2m r 



gk B Ti] 



(49) 



6 



Thus one obtains an infinite set of equations correspond- 
ing to the various power of ft 



iE ,o {0)aa ' 

A JP (n+l)aa' 



= o 



(50) 



/■ r N , \ (n)aa' 

E J -'^'PN e )/3/3 ' (n>l)-(51) 

00' 



As shown in Ref. [25|, in order to ensure that a solution 
can be found by recursion, one must discuss the solution 
of Eq. (|5l)l when calculating the diagonal elements 

To this end, 



in terms of the off-diagonal ones p^, " ■ 



using p^J 



J * aa ..f)>l3 and the 



Jaa,0f3' 

fact that J aa ,(3p = when a real basis is chosen, it is 
useful to re-write Eq. (|51|l in the form 



I ■ T N 1 \ l' 1 

(iL aa + n N )p y N 



{n)oiOL 



E 

0>0' 



91? ( T J n )00' 



(52) 



One has {—iL^ a — knV = iL^ a - The right hand side 
of this equation is expressed by means of the general- 
ized bracket in Eq. 1)26(1: and any general function 
f(H^) are constants of motion under the action of iL^ a . 
The phase space compressibility associated with the 
generalized bracket in the case of Nose dynamics is 




T • 



(53) 



where N is the number of classical momenta P in the 
Hamiltonian. Because of the presence of a non-zero phase 
space compressibility, integrals over phase space must be 
taken using the invariant measure |26| 



dM. = cxp(—w^)dRdPdrjdp n , 



(54) 



where u>§ = J dtn^ is the indefinite integral of the com- 



pressibility. To insure that a solution to Eq. (|52() exists 
one must invoke the theorem of Fredholm alternative, 
requiring that the right-hand side of Eq. I152H be orthog- 
onal to the null space of (iL^ a )^ [53 . The null-space of 
this operator consists of functions of the form 01 /(-^n)' 
where f(H^) can be any function of the adiabatic Hamil- 
tonian H?i. Thus the condition to be satisfied is 



[ dMj2 2 K ( J aa 

0>0' 



(n)00' 
00'PNe 



f(H$) = . (55) 



Apart from the integration on the additional Nose phase 
space variable there is no major difference with the proof 

given in Ref. (25|: 2K (,J a a,i30> P^f 13 ) and f(H§) are 
respectively an odd and an even function of P; this guar- 
antees the validity of Eq. (|55|l . 



Thus one can write the formal solution of Eq. I|52|l as 

(n)aa / ■ r N , \— 1 \ " on ( l i n )0 
Pke = ( lL aa + K N) > 2TZ J aa ,00'Piie 



0>0' 



and the formal solution of Eq. 115 1(1 for a ^ a' as 



(56) 



-(iL 



N 

act' 



\ (n)aa' 
KN)PNe 



Eaa ' 130' 



Jaa',Pf3'P Ne 



(57) 



Equations l|56|) and l|57|l allows one to calculate p^e to 
all orders in ft once p^"" is given. This order zero term 

is obtained by the solution of (iL^ a + ^rOp^"" = 0. All 
higher order terms are obtained by the action of E aa i, 
the imaginary unit i and J aa '/3p' (involving factors of 
d aa ' , P and derivatives with respect to P. Hence, one 
can conclude that functional dependence of p{J] m on the 
Nose variables rj and p n is preserved in higher order terms 

(n)aa f 
PNc ■ 

One can find a stationary solution to order ft by consid- 
ering the first two equations of the set given by Eqs. H5U|) 
and 



^N,P Ne 
#N,PNe 



= 



0), 



= -- (i?NA N /?N 



(0) 



~(°) A £rN 



(n = 1) . 

For the 0(ft°) term one can make the ansatz 



where Z is 



and obtain 



Z = ^2 J dM6(C-H§) 



(58) 



(59) 



(60) 



(61) 



PNe 



.P. ,„ 

-i-^d af) p' N 



(0)0 



. e -0(E a -E f3 ) 

Ep — E a 



2 



1 + e -0(E a -E ) 



(62) 



for the 0(ft) term. 

Equations H6U|) and l|62|l give the explicit form of the 
stationary solution of the Nose-Liouville equation up to 
order 0(ft). One can now prove that, when calculating 
averages of quantum-classical operators depending only 
on physical phase space variables, Q a (R,P), the canoni- 
cal form of the stationary density is obtained. It can be 
noted that it will suffice to prove this result for the 0(ft ) 
term since, as discussed before, the differences with the 
standard case are contained therein. 
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Indeed, when calculating 

(G a (R,P)) oc =J2 fdMg a (R,P) 

x5(C - H% - gk B T<n) . (63) 

Considering the delta function integral over Nose vari- 
ables, one has 



dp v dr) e- Nr] 6(C 



- H a - gk B Tr)) 

= const x exp[-(3(N/g)H^(R,P)] , 

(64) 



where it has been used the property S(f(s)) = 
[df /ds]j : } So S(s — so) (so is the zero of /(s)). Thus, at 
variance with what found in Ref. |16| . in order to re- 
cover the canonical distribution in the quantum-classical 
case, one must set g — N as it is done in the classi- 
cal case 0, 0- If the dynamics is ergodic and if one 
could integrate quantum-classical equations of motion for 
sufficiently long time, the phase space integral could be 
substituted by a time integral along the trajectory |28j . 
Ergodicity could be enforced by modern advanced sam- 
pling techniques |29j but long time stable integration of 
quantum-classical dynamics is still a challenge. 



For example, thermostated dynamics can be useful for 
preparing systems into desired initial conditions or for en- 
suring a good thermalization of the classical bath degrees 
of freedom providing a way to control the nonadiabatic 
character of the dynamics. On the conceptual side one 
could note that, historically, deterministic dynamics with 
thermodynamic constraints for purely classical systems 
has provided well defined algorithms to treat open sys- 
tems both in and out of equilibrium. Thus, the possibility 
to use the same tool in the case of quantum-classical sys- 
tems could disclose novel routes to the numerical study 
of open quantum systems. 

In conclusion, the non-Hamiltonian quantum formal- 
ism introduced in this paper sets a unified framework 
with the non-Hamiltonian classical algebra and, at the 
same time, discloses various routes for investigating gen- 
eralized quantum and quantum-classical systems. Such 
studies will be performed in the future. 
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VI. CONCLUSIONS AND PERSPECTIVES 



APPENDIX A: NHC AND NPT 



In this paper a generalized non-Hamiltonian form of 
quantum mechanics has been presented. This has been 
achieved through the introduction of a suitable non- 
Hamiltonian commutator which has been obtained by 
generalizing the symplectic structure of the standard 
quantum mechanical commutator. Therefore, it has been 
demonstrated that a single idea (i.e. generalizing the 
symplectic structure of the bracket while retaining its 
antisymmetric form) is able to describe in a unified way 
non-Hamiltonian theories both in classical and quantum 
mechanics. The non-Hamiltonian form of quantum me- 
chanics here presented provides a general mathematical 
structure which encompasses the ideas proposed by Wein- 
berg to introduce non-linear effects in quantum mechan- 
ics and whose physical content remains yet to be unveiled. 

For the sake of illustrating the possible use of non- 
Hamiltonian commutators, it has been shown that they 
subsume the quantum-classical bracket proposed by 
other authors. Moreover, their matrix structure has been 
used to define Nose dynamics on the classical degrees 
of freedom in quantum-classical systems. It has been 
also shown that the non-Hamiltonian quantum-classical 
bracket can be easily generalized to treat other thermo- 
dynamic constraints such as those provided by barostats 
or Nose-Hoover chains. The respective stationary den- 
sity matrices arc easily derived. The implementation of 
thermodynamic constraints for the classical degrees of 
freedom in quantum-classical systems could be consid- 
ered both as a practical and a conceptual improvement. 



The calculations of the previous sections show that the 
introduction of extended system dynamics on the classi- 
cal part of the system amounts to modify the operator 
of Eq. (|18|l by simply substituting the classical bracket 
operator given in Eq. I|16fl with the one suited to express 
the desired extended system dynamics [9L TTpj . 

Thus, in order to couple a Nose-Hoover chain to 
the classical coordinates, the classical phase space point 
is defined as 



(Al) 



where for simplicity one is considering a chain of just two 
thermostat coordinates r\\ , 772 and momenta p m , p ri2 . 



H 



NHC 



P 2 



2m 2M 2 

$(<?, R) + gk B T m + gk B T m , (A2) 



where m. ni and vn m are the inertial parameters of the 
thermostat variables. As shown in Ref. 0,^3 j one can 
define an antisymmetric matrix 



NHC 



" 








1 





" 














1 




















1 


-1 











-p 








-1 





p 





-Pvi 








-1 





Pvi 






(A3) 
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The matrix j£? NHC determines the operator A NHC which 
in turn provides the the non-Hamiltonian bracket accord- 
ing to Eq. The Nose-Hoover chain classical equa- 
tions of motion in phase space are then given by 



X = -XA NHC i? NHC . 



(A4) 



Quantum-classical dynamics is then introduced using the 
matrix super-operator 



V 



NHC 







1 + A A NHC 



(1 + ^ANHC) 







(A5) 



As previously shown by means of the latter the quantum- 
classical equations of motion are then given by 



Q = i [ NHC 

dt h L 



NHC 



X 



(A6) 



The equations of motion can be represented using the 
adiabatic basis obtaining the Liouville super-operator 



t nhc 

lJ ~aa',/3(3' 



(A7) 



where 



The equations of motion for the classical coordinates 
are @ 



with 



R 

fin 
V 

P 

fin 
fiv 



mi 
Pv 
m. n 
Pv 
my 
<93> 



R 



Pv 



3Vmv 



p Pv 

dR 3Vm v 

1 P 2 P\ 



P 



Pi! 



= £-+ — ~9k B T 
Vfli my 

2—1 

T-i PV 

= F v - pv — 
m v 



F v = — 



N 



2-i M dR 



(A12) 
(A13) 
(A14) 
(A15) 

(A16) 
(A17) 

(A18) 



The antisymmetric matrix to define the Operator A NPT 
in Eq. (jT1)|) is then 



iL 



NHC 



P_d_ 

MdR + 

Pm 



1 



(F a + F a ) 



dP 



2 

£■ 

k=l 

Pm 
m 







d 







dp 



1)2 U 1 J V1 



(A8) 



with F V2 = p^ li /m r]1 — gksT. The proof of the existence 
of stationary density matrix in the case of Nose-Hoover 
chains follows the same logic of the simple Nose-Hoover 
case. In the adiabatic basis the density matrix stationary 
up to order bar has the same form as given in Eqs. i|(j(J|) 
and H62|) . One has just to replace Eq. (|60|l for the order 
zero term with 



(0)q/3 
rNHCe 



1 -P 



- +E a (R) 2 k= i ( T^h- +gk B T m 



(A9) 



with obvious definition of Z. 

For the case of constant pressure and temperature dy- 
namics, the equations of motion treated in Ref j^lj are 
here considered. This time, the extended phase space 
point is 



B 



NPT 







-1 



R 
3V 

2 NPT 








-1 

-1 



1 





p 

p 

3V 





1 



-p 



-Pv 



R 

3 V 



1 

p 

3V 

Pv 




(A19) 



The matrix is the same as that given in Ref. 

but this time the order of coordinates and momenta in 
the classical extended phase space point definition of 
Eq. (| A10|) ensures that A NPT in Eq. (jl6|l makes the equa- 
tions of motion 



X 



-XA 



NPT^NPT 



(A20) 



exactly equivalent to Eqs. (|A12IA17|) . Then one can 
define the matrix operator T) and write down the 
quantum-classical equations of motion. In the adiabatic 
basis the equations are written by means of the Liouville 
operator 



r NPT 



(iu) aa > + iL aa , )6 af 35 a >f3> 



(A21) 



with 



X = (R, V ,V,P,p v ,p v ) , (A10) 
and the Hamiltonian quantum-classical operator is 



H 



NPT 



f_ P*_ 

2m + 2M 



Pv 



2m r 



i-jjj 2mv 
&(q, R) + gk B Tr, + P cxt V . (All) 



■ r NPT _ ■ r NH 



Pv -R 9 



Pv p d 



Wm v dR 3Vm v dP 



Pv d p v d 

mv dv 



m n dP\ 



V 



(A22) 



The stationary density matrix is derived as usual and in 
the adiabatic basis is expressed again in the form given 
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by Eqs. and © with the O(h ) term given by 

(0)a/3 _ 1 -p\&+EaW+^+gk B Tr,+ £jfc+P. m tV 



(A23) 
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